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\—i ■ Abstract. We show that the Euchdean Wasserstein distance is contractive for inelastic 

^~H| homogeneous Boltzmann kinetic equations in the Maxwelhan approximation and its asso- 

ciated Kac-like caricature. This property is as a generaHzation of the Tanaka theorem to 
inelastic interactions. Even in the elastic classical Boltzmann equation, we give a simpler 
proof of the Tanaka theorem than the ones in J25L I77| . Consequences are drawn on the 
asymptotic behavior of solutions in terms only of the Euclidean Wasserstein distance. 

<^ ■ 1. Introduction 

■5 
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This work is devoted to contraction and asymptotic properties of the homogeneous 
Boltzmann-type equations for inelastic interactions in the Maxwelhan approximation in- 
troduced in j3] and further analyzed in jISl El El 13 IH CHI 12 E] • We are basically concerned 
with the Boltzmann equation 

c^' % = BVo{mQ{fj) (1-1) 

m '■ 

en ' considered in [^ and its variants. Here, f{t,v) is the density for the velocity t; G M'^ 

^ ■ distribution of the molecules at time t, and Q{f,f) is the inelastic Boltzmann collision 

\^ ! operator defined by 

;:^: (^^Q(fj)) = Lf f [ f{v)f{w)L{v')-v{v)]dadvdw (1.2) 

^ ■ 47r J^3 J^3 Js-2 L J 



S 



for any test function Lp, where 



K.^ /!/ xl^e, ,1+e, , 

•^ ■ V = -[v + w) A [v — w) ^ : — \v — w\a 

. r-j . 2 4 4 

j_| ■ is the postcollisional velocity, a G 5^, f,w G M^ and < e < 1 is the constant restitution 

^ coefficient. Equation (jl.lj) preserves mass and momentum, but makes the kinetic energy 

(or temperature) 

Oifit)) = \[ V- [ vf{t,v)dv'f{t,v)dv 

•J Jm.3 Jm.3 

decrease towards 0. In particular, solutions to ()1.H) tend to the Dirac mass at the mean 
velocity of the particles jH] . We refer to jiSl El I2H1 for the discussion about the relation of 
this model to the inelastic hard-sphere Boltzmann equation and different ways of writing 
the operator. Let us just point out that the factor B \JO{f{t)) in front of the operator in 
()1.1|) is chosen for having the same temperature decay law as its hard-sphere counterpart 
[^ known as the Haff 's law. 

1 



2 FRANgOIS BOLLEY AND JOSE A. CARRILLO 

The convergence towards the monokinetic distribution has been made more precise in 
[71 ini 121 by means of homogeneous coohng states. They are self-similar solutions of the ho- 
mogeneous Boltzmann equation (jl.lll describing the long-time asymptotics and presenting 
power-like tail behavior whose relevance was previously discussed in the physics literature 
[16, 17J. 

To avoid the collapse of the solution to the Dirac mass, the authors in ^^ suggested the 
introduction of a stochastic thermostat which, at the kinetic level, is modelled by a linear 
diffusion term in velocity. In this framework, the density / in the velocity space obeys 

^ = Bye{fW)QifJ) + Ae^{f{t))Aj with o<p<^- (1.3) 

Existence and uniqueness for given mean velocity of a steady state to ()1.3|1 have been shown 
in [iniinilll- The convergence of solutions towards this steady state in all Sobolev norms 
has also been investigated and quantified by means of Fourier-based distances between 
probability measures [1^ 

Fourier techniques are a good toolbox and have been extremely fruitful for studying 
Maxwellian models in kinetic theory since Bobylev observed j2l E] that such equations 
have closed forms in Fourier variables. Fourier distances are not only suitable technical 
tools to study the long-time asymptotics of models p.lj) and p.Hj) . but also they represent 
the first Liapunov functionals known for inelastic Boltzmann-type equations Pl E]- In 
the case of the classical elastic Boltzmann equation for Maxwellian molecules, there is 
another known Liapunov functional, namely, the Tanaka functional |2S], apart from the 
if -functional for which no counterpart is known in inelastic models. 

The Tanaka functional is the Euclidean (or quadratic) Wasserstein distance between 
measures in the modern jargon of optimal mass transport theory. It is defined on the set 
V2{^^) of Borel probability measures on M.^ with finite second moment or kinetic energy 
as 

W2{f,g) = mU [[ \v-w\^d7r{v,w)\ = ini {E[\V -W\^]y^^ 

where vr runs over the set of joint probability measures on M^ x M^ with marginals / 
and g and (V, W) are all possible couples of random variables with / and g as respective 
laws. This functional was proven by Tanaka [25j to be non-increasing for the flow of the 
homogeneous Boltzmann equation in the Maxwellian case. In fact, the Tanaka functional 
and Fourier-based distances are related to each other [THl [12112^1, and were used to study the 
trend to equilibrium for Maxwellian gases. On the other hand, related simplified granular 
models ^1] have been shown to be strict contractions for the Wasserstein distance W2- 

With this situation, a natural question arose as an open problem in Remark 3.3] 
and [2H1 Section 2.8]: is the Euclidean Wasserstein distance a contraction for the flow of 
inelastic Maxwell models? The main results of this work answer this question affirmatively. 
Moreover, we shall not need to introduce Bobylev's Fourier representation of the inelastic 
Maxwell models working only in the physical space. 

We shall show in the next section the key idea behind the proof of all results concerning 
contractions in W2 distance for inelastic Maxwell models, namely, the gain part Q~^{f,f) 
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of the collision operator verifies 



W2{Q^{f,f),Q^{g,g))<^^^W,{f,g) 

for any f,g in V2{^'^) with equal mean velocity and any restitution coefficient < e < 1. 
Based on this property, we shall derive contraction and asymptotic properties both for 
(jl.lll and p.Hj) in Subsections Kill and IH.2I On one hand, we shall prove that the flow for 
the diffusive equation ()1.H|1 is a strict contraction for W2, while for the scaled equation 
associated to (jl.lj) we shall show that solutions converge in W2 to a corresponding ho- 
mogeneous cooling state, without rate but only assuming that initial data have bounded 
second moment. This improves the Ernst-Brito conjecture fHl ^1 HI HI 12] since it shows 
that the basin of attraction of the homogenous cooling state is larger -we avoid the typical 
assumption of bounded moments of order 2 + 5- if we do not ask for a rate. 

Moreover, a generalization for non constant cross sections including Tanaka's theorem 
as a particular case will be proven in Section EJ Finally, we shall also show this generic 
property for the inelastic Kac model introduced in [221 ^^ ^ dissipative version of Kac's 
caricature of Maxwellian gases [T^ I2Uj. 

2. Contraction in W2 of the gain operator 

We start by summarizing the main properties of the Euclidean Wasserstein distance W2 
that we shall make use of in the rest, refering to ^U |2Zj for the proofs. 

Proposition 1. The space (7^2 (I^^), 1^2) 'is a complete metric space. Moreover, the fol- 
lowing properties of the distance W2 hold: 

i) Convergence of measures: Given {/„}«>! O'nd f in V2{^^), the following three 
assertions are equivalent: 

a) W2{fn, f) tends to as n goes to infinity. 

b) /„ tends to f weakly-* as measures as n goes to infinity and 



sup 



■P / \'v\^ fn{v) dv -^ Q as i? ^ +00. 

:1 J\v\>R 

c) /„ tends to f weakly-* as measures and 

\v\^fn{v)dv^ I \v\^f{v)dv as n ^ +00. 



iii) Relation to Temperature: /// belongs to ^2(1^^) o-nd 5a is the Dirac mass at a 
in M.^ , then 

W^if,5a)= [ \v-a\'df{v). 
iii) Scaling: Given f in V2{^^) and 6 > 0, let us define 



Se[f] = ey'fie'/\) 
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for absolutely continuous measures with respect to Lehesgue measure or its corres- 
ponding definition by duality for general measures; then for any f and g in V2{^^), 
we have 

W2{Se[f],Se[g])=9-'/'W2{f,g). 
iv) Convexity: Given fi, f2, gi and g2 in V2(^^) and a in [0, 1], then 

W^iaf, + (1 - a)f2, ag^ + (1 - a)g2) < aW^{f,, g,) + (1 - a)W^{f2, ^2). 

As a simple consequence, given f,g and h in V2{^^), then 

W2{h*f,h*g)<W2{f,g) 

where * stands for the convolution in M.^. 

Here the convolution of tlie two measures h and / is defined by duality by 



{^,h*f)= ip{x + y)dh{x)df{y) 

for any test function y? on M'^. If / is a Borel probability measure on M^ we shall let 

</>= / vdf{v)= / vf{v)dv 



denote its mean velocity. We shall use the same notation for densities and measures 
expecting that the reader will not get confused. 

Let us write the collision operator Q given in ()1.3|) as 

Q{fJ) = Q^{fJ)-f (2.1) 

where Q~^{f, f) is defined by 

{V.Q^UJ)) = ^I I f f{v)f{w)^{v')dadvdw (2.2) 

47r J^i Jr3 Js2 

for any test function (f, where we recall that 

, 1 , , 1 - e , , 1 + e, , 

V = -[V + w) -\ : — [V — w) -\ : — \v — w\a. 

2 4 4 

In this section we derive a contraction property in W2 distance of the gain operator 

Q^. For that purpose, let us note that the previous definition of the gain operator can be 

regarded as follows: given a probability measure / on R^, the probability measure Q^{f, f) 

is defined by 

(¥',Q+(/,/))= / I f{v)f{w){ip,Yi,^^)dvdw 



where Iiv,w is the uniform probability distribution on the sphere S^^w with center 

1, , 1-e 

-[v + w) + — 

and radius 



(^v,w = o(^ + w)^ ^— (t; - w) 



1 + e, , 

rv,w = -^—\v-ui\. 
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Figure 1 . Geometry of inelastic collisions 

In probabilistic terms, the gain operator is defined as an expectation: 

Q+{f,f)=E[Uv,w] 

where V and W are independent random variables with law /. 
Then the convexity of W2 in Proposition ^ implies 

<E[W^{nv,w,'nx,Y)] (2.3) 

where X and Y are independent random variables with law g. This observation leads us to 
consider the W2 distance between uniform distributions on spheres. To this aim, we have 
the following general lemma: 

Lemma 2. The squared Wasserstein distance W2 between the uniform distributions on 
the sphere with center O and radius r and the sphere with center O' and radius r' in M? is 
bounded by \0' — Op + (r' — r)^. 

Proof.- We define a map T : M^ — > M.^ transporting the sphere of center O and radius 
r > onto the sphere with center O' and radius r' > r in the following way: 

• If r = r', then we just let T be the translation map with vector O' — O, i.e., 
T{v) =v + 0'-0. 

• If O = O', then we just let T be the dilation with factor — centered at O, i.e., 
T{v) = ^^v. 

• If r 7^ r', then we consider the only point fi G M^ verifying that 



that is, 



Q = + 






-iO'-O). 
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Then we let T be the dilation with factor — centered at Q, that is, we let T{v) = 
Q + —{v — Q). Such a construction of the point Q and the map T is sketched in 
Figure |21 in the case of non interior spheres. 




T(v) 



Figure 2. Sketch of the computation of the Euclidean cost of transporting 
spheres to spheres. Transport lines are just rays from the point fl. 



Let h(o,r and Uo'y denote the uniform distributions on the corresponding spheres. Then 
the transport plan it given by 



ri{v,w) d'7i{v,w) 



r]{v,T{v))(HAo,r{v) 



for all test functions ri{v,w) has Uo,r and Uo'y as marginals by construction of T. Using 
this transference plan in the definition of the Euclidean Wasserstein distance, we finally 
conclude 



W^{Uo,r,Ko',r') < 



\v-T(v)\''dUo,rM 



I _ \ 2 



k' — fiP dUoriv) 



that can be computed explicitly, giving 

W^{Uo,r,l^o'y) < \0' - Op + (r' - r^ 

and finishing the proof. D 

This lemma, using the notation a = v — x and b = w — y, ioi fixed values v, w, x, y in 
M^, implies that 



W^{U^^^,U^^y) < \a 



v,w ^x,y\ 



< 



+ |r 
1 + e 



v,w 'x,y 



3-e 

-^ ^+ A 

4 4 

5 - 2 e + e^ 



+ 



l + e\2 



aP + 



4 



|a-6p 



1 -p2 



here a ■ b denotes the scalar product between a and b in M^ and the bound in 



I ' v,w ' x,y\ 



< 



1 + e 

4 
l + e\2 



\v — w\ — \x — y\ 



J \{v ~ w) - {x - y)\ 



l + e\2 



) l«-&l' 
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follows from the Cauchy-Schwarz inequality 

{v — w) ■ {x — y) < \v — w\ \x — y\. (2.4) 

Therefore, by ()2.3|) . 

W^{Q^f,f),Q^{g,g)) <^-^^^E[\V - X\'] + ^^^E[\W -Y\'] 

Let moreover (V, X) and {W, Y) be two independent optimal couples in the sense that 

W^if, g) = E [\V - X|2] =E[\W- Fp] . 
Then 

E[{V -X)-{W -Y)] =E[{V-X)] ■ E[{W -Y)] = \ </> - <g> |^ 

by independence. Collecting all terms leads to the following key estimate and contraction 
property: 

Proposition 3. /// and g belong to V2{^^), then 

WiiQ^ifJ),Q^ig,g))<^^Wi{f,g) + ^-^\<f>-<g>\' 

for any restitution coefficient < e < 1. As a consequence, given f and g in V2{^^) with 
equal mean velocity, then 



W2iQ'-{fJ),Q^{g,g)) < ^j^-^W^if.g). 

The case of equality is addressed in the following statement: 

Proposition 4. Let f and g belong to V2{^^) with equal mean velocity and temperature, 
where g is absolutely continuous with respect to Lebesgue measure with positive density. If 



W2iQHfJ),QH9,9)) = ^^^W2if,g). 
for some restitution coefficient < e < 1, then f = g. 

Proof.- It is necessary that the equality holds at each step of the arguments in Proposi- 
tion |21 In particular, (j2.4p holds as an equality, that is, 

V-W _ X-Y 

\V-W\ ~ \X-Y\ 

almost surely in the above notation. Then, since g is absolutely continuous with respect 
to Lebesgue measure with positive density, one can proceed as in |25[ Lemma 9.1] to show 
that f = g. We sketch the proof for the sake of the reader. Since g is absolutely continuous 
with respect to Lebesgue measure, there exists [2Z] a Borel map -u : R'^ — i> R^ such that 
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/ be the image measure of g by u, and in probabilistic terms V = u{X) and W = u{Y) 
almost surely. Hence 

u{x) — u{y) X — y 



\u{x)-u{y)\ \x-y\ 



(2.5) 



almost everywhere for Lebesgue measure since X and Y are independent and since their 
law g has positive density. We leave the reader to check |27| Exercise 7.25] that this implies 
the existence of constants ui and 002 such that u{x) = uji + UJ2X. First of all cug = 1 since 
/ and g have same temperature. Then identity (j2.5p forces a;2 = 1, implying cui = since 
</>=<(?>, and finally / = 5^. D 



3. Contractive Estimates for the Inelastic Maxwell Model 

In this section, we shall derive contractive estimates in the Euclidean Wasserstein dis- 
tance for solutions to the inelastic Maxwell models both in the non-diffusive and the dif- 
fusive cases. 

3.1. The non-diffusive case. We are first concerned with solutions /(t) to the Boltzmann 
equation p.l|) with < e < 1. After time scaling defined by 

Q 

with E = -, as in [?, we get a function denoted again /(r) for simplicity, solution to 

1 — e^ 

^ = EQ{f,f). (3.1) 

Theorem 5. If fi and /2 are two solutions to ()3.Hl with respective initial data /° and /^ 
%n V2iM^), then 

W^|(/i(r),/2(r)) < e-2- iy|(/°, /°) + (1 -e'^^ |</°> - <f^>\' (3.2) 

for all T > 0. 

Proof.- Decomposition (j2.ip of the collision operator Q as 

Q{fJ) = Q^{fJ)-f 

allows us to represent the solutions to ()3.1|) by Duhamel's formula as 

/.(r)=e-^-/° + E re-^(--)Q+(/,(.), /,(.)) rfs, ^ = 1,2. 
Jo 
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Then the convexity of the squared Wasserstein distance in Proposition ^ and Proposition 
121 imply 

W|(/i(r),/2(r)) 

< e-^- Wiifl f,) + E /V^(--) Wi (Q^Ms), Ms)), Q+(/2(5), f2{s))) ds 

Jo 

< e-^^ W^K/?, /°) + E £e-^(^-) (^^ W^{Ms),Ms)) + X^ ds; 



here 



X = ^ I <flis)> - <f2is)> \' 



does not depend on time since the mean velocity is preserved by equation ()3.ip . In other 
words, the function y^r) = e^'^ W^2^(/i(r), f^ij)) satisfies the inequality 

y{r) < y(0) + E j^^ {^-^ y{s) + X e^^^ ds 

and then 

l/(r)<i/(0)e^^^ + -^(e^--e^^-) 

1-7 

by Gronwall's lemma with 7 = (3 + e^)/4. This concludes the argument since (1— 7) E = 2. 
U 

Remark 6. 

(1) Without further assumptions on the initial data /° and /^ , this result is optimal in 
the following sense. If /^ is chosen as the Dirac mass at the mean velocity of f^, 
then inequality ()3.2|) is actually an equality for all r; indeed 

W^{Mt),Mt))= f \v-<Mr)>\'MT,v)dv = 3e{Mr)) 

JIR3 

= 3e-2-^(/{')=e-2-W^|(/°,/°) 

JO 

since -^ = — 2 6* by equation ()3.1|) . 

(2) In terms of the original time variable t in (II. Ij) . if /{* and /j are two initial data 
with the same initial temperature ^O) then the temperatures of the corresponding 
solutions /^ and p to ()1.1|) follow the law 



d9 1-e 



,2 



,3 



dt 4 

and hence are both equal to 



B92 (3.3) 



3-1/2 



1 2 \ -2 
1 — e^ 



^(t) = 00 '' + — — 5t 
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Then estimate ()3.2|) reads as 



for all t > 0. 



w^i(/r,/2°) + (i-^)i</?>-</2°> 



The convergence of the solutions to (ll.lj) towards the Dirac measure at their mean velocity 
has been made precise in P[^ by the introduction of self-similar variables and homogeneous 
cooling states. There the authors prove that the rescaled solutions g defined by 

g{T,v)=9'/\f{T))fir,9'/\f{T))v) (3.4) 

satisfy the strict contraction property 

4+.(^i(r),^2(r)) < e-^(^)-rf2+.(^?,^2°), Cie) > 

for initial data g^ and g2 in V2{^^) with equal mean velocity and pressure tensor, where 
e > and d2+e is a Fourier-based distance between probability measures. Moreover, for 
e = one has C{e) = giving a non-strict contraction in c?2 distance. In fact, by the 
scaling property in Proposition [H (j3.2|) reads as 

W2igiiT),g2iT))<W2iglg'2) (3.5) 

in the scaled variables. This is consistent with the fact that the distances d2 and W2 are 
"of the same order" [THl I2ni 12] up to moment bounds. 

A measure g{T, v) defined by (j3.4j) from a solution f{v, r) to (|3.1|) with initial zero mean 
velocity has zero mean velocity and unit kinetic energy for all r, and is solution to 

^ + V-{gv) = EQ{g,g). (3.6) 

Moreover it is proven in [HI |2] that (|3.6|) has a unique stationary solution g^o with zero 
mean velocity and unit kinetic energy; all measure solutions g{T,v) to (j3.6|) with zero 
mean velocity, unit kinetic energy and bounded moment of order 2 + 6 converge to this 
stationary state goo as r goes to infinity in the d2 sense, that is, in the W2 sense since d2 
and W2 metrize the same topology on probability measures [22] up to moment conditions. 
Moreover the convergence has exponential rate in the ^2 sense, and in the W2 sense if the 
initial datum has finite fourth order moment. In turn this ensures existence and uniqueness 
of homogeneous cooling states to (jl.l|) for given mean velocity and kinetic energy, and 
algebraic convergence of the solutions f{t) towards them in the original variables. 

We conclude this section by proving this convergence result using only the W2 distance, 
and without assuming that the initial data has more than two finite moments. This in turn 
shows that the Euclidean Wasserstein distance W2 between solutions of ()3.6|) converges to 
zero as t goes to infinity, improving over ()3.5|) that does not a priori yield any information 
on the long-time behavior of the solutions g. As a drawback, this argument does not 
provide any rate of convergence as does the Fourier-based argument in j2j . 
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Theorem 7. Let g^ and g2 be two Borel probability measures on M'^ with zero mean velocity 
and unit kinetic energy, and let (71 (r) and (72 (t) be the solutions to (j3.6p with respective 
initial data g^ and g2- Then the map t i-^ VF2(5'i(t),5'2(t)) is non-increasing and tends to 
as T goes to infinity. 

Proof.- It is based on the argument in ^7\ to Tanaka's theorem. The first statement is 
a simple consequence of ()3.5|) . Then we turn to the second part of the theorem which by 
triangular inequality for the W2 distance is enough to prove when g2, and hence (72 (t), is 
the unique stationary state goo to ()3.6p with zero mean velocity and unit kinetic energy. 

Step 1.- Let us first assume that the fourth moment of the initial datum is bounded, 
i.e., 

/ l'^l'^5'i('^) dv < 00. 

Then Proposition ^J in the appendix ensures that 

sup / Ivl"^ gi{T,v) dv < 00, 



T>0 

so that 

sup / \v\'^ gi{T,v) dv 

T>0 J\v\>R 

tends to as i? goes to infinity. Prohorov's compactness theorem and Proposition pimply 
the existence of a sequence r^ ^ 00 as A; ^ 00 and a probability measure /i° on R^ with 
zero mean velocity and unit kinetic energy such that W2{gi{Tk) , fJ^^) —>■ as k -^ 00. We 
want to prove that /i° = goo- 

Without loss of generality, we can assume that the diverging time sequence satisfies 
'i'fc + 1 < Tfc+i for all k. Now, since goo is a stationary solution, it follows from the first part 
of the theorem that 

W2{gi{Tk+i),goo) < W2igiirk + l),goo) < W2{gi{Tk),goo)- (3.7) 

On one hand, both W2{gi{Tk),goo) and W^2(fi'i(T"fc+i),fl'oo) tend to W2{fi°,goo) as k goes to 
infinity by triangular inequality. Then, if /i(r) denotes the solution to (|3.6|) with initial 
datum /x°, the first point again ensures that 

1^2(^1 (Tfc + l), Ml)) <W^2(^l(rfc),/i°) 

which tends to 0. Hence W2{gi{Tk + l),goo) tends to W2{fi{l),goo) by triangular inequality, 
and finally 

W2i^Iil),goo) = W2i^^',goo) 
by passing to the limit in k in ()3.7|) . By the non- increasing character of W2 along the flow, 
we deduce that 

W2{l^il),goo) = W2{Kr),9oo) = W2{lX^goo) 

for all T e [0,1]. 
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Consequently /i(r) and g^o are two solutions to ()3.6|) with zero mean velocity and unit 
temperature, whose W2 distance is constant on the time interval [0,1]. This is possible 
only if equality holds at each step in the proof of Theorem El in the original space variables; 
in particular 



W2(Q+(/iW,/i(r)),g+(^oo,^oo)) = ^J^-^W2{Kr),9oo) 

for all r, and especially for r = 0. But fi^ and Qoo have same mean velocity and temperature, 
and, according to [3 Theorem 5.3], f^oo is absolutely continuous with respect to Lebesgue 
measure, with positive density. Hence Proposition E] ensures that /i° = Qoo- 

In particular W2{gi{Tk),goo) -^ as /c — > 00, and then W2(5'i(t),5'oo) — > as r — i> 00 
since it is a non increasing function. 

Step 2.- Let us now remove the hypothesis on the boundedness of the initial fourth 
order moment. Let {g^"')n be a sequence in ^2(1^^) with zero mean velocity, unit kinetic 
energy, finite fourth order moment and converging to Qi in the weak sense of probability 
measures; in particular it converges to g^ in the W2 distance sense by Proposition [TJ Such 
a g^"" can be obtained by successive truncation of g^ to a ball of radius n in M'^, translation 
to keep the mean property, and dilation centered at to keep the kinetic energy equal to 
1. 

Then, if (/"(r) is the solution to (j3.6p with initial datum g^"-, the triangular inequality 
for W2 and (j3.5p ensure that 

W2{gi{T),g^) <W2{gi{T),g^{T)) + W2{g''{T),g^) 

<W2{g'^,g'n+W2{g''{r),goo). 

Given e > 0, the first term in the right hand side is bounded by e for some n large enough, 
and for this now fixed n, the second term is bounded by e for all r larger than some 
constant by the first step. This ensures that W2{gi{T),goo) tends to as r goes to infinity. 
D 

3.2. The diffusive case. We now turn to the diffusive version ()1.3j) of p.ll) . Again by 
the change of time 

Q 

with E = we are brought to studying the equation 



1-e' 



^ = EQ{fJ) + &\f{T))AJ {3.i 



where 



^"^-imffr)]]"-'^' 
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As in the nonviscous case of ()3.H) we shall prove 

Theorem 8. If fi and /2 are two solutions to ()3.8|) for the respective initial data /° and 
/2 in V2{^^) with same kinetic energy, then 

W^I(/iM,/2(r)) < e-^-iyK/?,/") + (1 -e-2-) |</°> - <f',>\' (3.9) 

for all T > 0. 

Proof.- We again start by giving a Duhamel's representation of the solutions. To this 
aim we write ()3.8)1 as 

EF-Ef + Q\fiT))Af 



dr 
where F = Q^{f, /), that is, 



^+[E + \k\'Q\f{T)))f = EF. 
Here, we are using the convention 

fx{k) = / e'"'''' diJ,{x) 
for the Fourier transform of the measure /x on M'^. Hence the solutions satisfy 

/(r, k) = e-^- f%k) e-^(^'-) 1^1' + E T e'^^^-'^ F(s, k) e-(2(/--)-^(/-^))l'=P ds 

Jo 

where S(/, r) = / Q'^{f{s)) ds, and thus 
Jo 

f{T,v) = e-^- if * r2^if,r)){v) + E /V^^--^) {F{S) * T2i^if,r)-m,s))){v) ds 

Jo 

:= e-^" /(r, v) + E f e-^^^-'^ F{t, s, v) ds. 
Jo 

Here 

r (v) = i e-l^l'/2a 

is the centered Maxwellian with temperature a/3 > 0. Moreover /i and /2 have same 
temperature at all times, so that S(/i,r) = E(/2,r). Then the convexity of the squared 
Wasserstein distance and its non-increasing character by convolution with a given measure, 
see Proposition n imply that 

W^I(/iM,/2(r))<e-^-iy|(A(r),/2(r)) + E/W-)iy2'(A(r,s),^2(r,.))rf5 

Jo 

< e-^- W^ifl /°) + E r e-^(-~^) W^{F,{s), F^^s)) ds. 

Jo 

In other words the squared distance lV2^(/i(r), /2(r)) satisfies the same bound as in the 
nonviscous case of Theorem and we can conclude analogously. D 
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Remark 9. 

(1) As pointed out to us by C. Villani the result can also be obtained by a splitting 
argument between the collision term and the diffusion term. 

(2) As proven in pQ, the temperature 0{f{t)) of the solution / in the original time 
variable t converges towards 

„ , 8 A \ 3^2^ 



.S(l-e2) 
as t goes to infinity, and satisfies 9{f(t)) > min(6'(/(0)), 6'oo). In particular 



r = §I^VWW)ds>^t 



iiCi = B min(^(/(0)), ^oo)^ • Writing ()H.9|1 in the original variable t for initial data 
with equal mean velocity and temperature, we recover the contraction property 

that coincides with (3.1) in Ij for the Fourier-based d2 distance exactly with the 
same rate. For p = 1 one can exactly compute r and also recover (3.2) in PP but 
for the distance W2- 
(3) The existence of unique diffusive equilibria for each given value of the initial mean 
velocity can be obtained from this contraction property of the W2 distance analo- 
gously to the arguments done in [T] with the Fourier-based distance ^2. 

4. General cross section 

In this section, we consider the more general case of a variable collision cross section 
when the gain term Q"*" is defined by 

{'P,Q^{f,f)) = ^ [ [ [ f{v)f{w)viv')b(^^^.cr)dadvdw 
4vr J]R3 Jk3 Js2 \\v~w\ J 

where again the post-collisional velocity v' is given by 

/ 1/ X l~e, , 1 + e, , 

V =-[v + w)-\ [V — w) -\ \v~w\a 

2 4 4 

and the cross section b satisfies the normalized cut-off assumption 

/ b{k-a)da= / 6(cos^) sin^t/^t/0 = 1 (4.1) 

^52 Jo Jo 

for any k in S^. Then we shall prove the following extension of Proposition |2] for non 
constant cross sections b: 

Theorem 10. /// and g inV2(^^) have equal mean velocity, then 
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Before going onto the proof, we draw the main consequence. Let /i = fi{T,v) and 
/2 = f2{'Ti'v) be two solutions to the Boltzmann equation 

^ = QU\f) = Q^U\f)-f 

with respective initial data /{* and /g in 7^2 (IR^), where Q^ is defined as above. Then, as 
in Section im Duhamel's representation formula 



/(r) = e-^ f + r e-(^-) g+(/(.), /(.)) ds 
Jo 



of the solutions and the convexity of W2 ensure the contraction property 

W^2(/i(r), /2(r)) < e-(^-^^)-/2 W^if,, f,) (4.2) 

for all r, where 

3 + e^ 1 - e^ r 
7fe = — \ — Tx I 6(cos 6) cos 6 sin ^ d^ 

4 ^ Jo 

is bounded by 1 by ()4.1|1 . 

In the elastic case when e = 1, 7;, = 1, one recovers Tanaka's non-strict contraction 
result |2^ for the solutions to the homogeneous elastic Boltzmann equation for Maxwellian 
molecules, at least under the cut-off assumption, but with a somehow simpler argument 
than those given in ^^ and 



Proof.- By definition 

(^,Q+(/,/)) = 2^ r [ [ \ [\{v')^}fiv)fiw)dvdwbicose)smede 

Jo Jr^ Jr3 wo ^^ -• 

= 27r / E[{(p,Uv,w,e)]b{cose) sinOde 
Jo 

where V and W are independent random variables distributed according to / and, given 

v,w in R^, Uv,wfi is the uniform probability measure on the circle Cy^w,e with center 

1/ N /1-e 1 + e A, , 

Cv,w,e = -{v + w) + (^— ^ — cos6'j(t;-w;), 

radius 

and axis 

, V — w 

k = :■ 



\v — w\ 
Let also g he a. Borel probability measure on M^ and X, Y be independent random 
variables with law g. Then, by the normalization assumption ()4.1|) . the convexity of the 
squared Wasserstein distance with respect to both arguments ensures that 

W^{Q+{f,f),Q+{g,g)) <27t [ E [W^{Uv,w,e,Kx,Y,e)] bicosO) sm9d9. (4.3) 

^0 
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Figure 3. 

We now let v,w,x,y and 9 be fixed in M.^ and [0,7r] respectively, and give an upper 
bound to W2(U^^w,eMx,y,e)- This consists in estimating the transport cost of a circle in M^ 
onto another one, for which we have the following general bound: 

Lemma 11. ^7\ The squared Wasserstein distance between the uniform distributions on 
the circles with centers c and c' , radii r and r' and axes k and k' is bounded by 

\c - c'p + r^ + r'^ - rr'(l + \k ■ k'\). 

Hence, using the notations a = v — x,b = w — y,d = v — w and b = x — y in our case we 
get 



W2(Uv,w,eMx,y,e) < 



< 



^-e l+e \ 1+e 2 

1 cos u ]a -\ ( 1 — cos a ) 

4 4/' 

sin^^ 
3-e 1+e 



l + e\2 . 2 



lap 



+ 2 



4 4 

l + e\2 



\b\'-Mb\[l+[--- 
2 /l + e\2 



a 



H '^- cos 6] + ( ^—^ ] sin^ 6 

3-e 



a\ |6|yJ 
2 



\a 



4 
l + e\2 



+ COS0 (1-cos^) -2sin2^ 



4 



1 + e 
[{l-cosef + sm^e]\b\ 



a ■ b 



(4.4) 



where we have used the bound 

|a|2 + |6|2 - \a\\b\ - a ■ b < lap + |6|2 - 2d ■ b = \d - b\'' = \a - ftp. 

Assume now that {V, X) and {W, Y) are two independent couples of random variables, 
optimal in the sense that 



Note that 



W^{f,g) = E[\V- X\'] =E[\W- Y\'] . 
E [{V -X)-{W -Y)]=E [{V -X)]-E [{W -Y)]=0 
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since (V, X) and {W, Y) are independent and since / and g have same mean velocity. Then 
from (jO|) : 

E[W^{Uv,w,e,Ux,Y,0)] <l{0)W^{f,g) 
where 



7W 



3 % 1±^ cos^)' + (^)'[(1 - cos^)2 + 2sin^ 



4 



4 



3 + e2 l-e^ 

H ^ — cos u. 



4 4 

One concludes the argument after averaging over 6 as in ()4.3|) and taking ()4.ip into account. 
D 



5. Inelastic Kac Model 

In this last section we consider a simple one-dimensional model introduced in [SB] which 
can be seen as a dissipative version of the Kac caricature of a Maxwellian gas f^ 120] • 
Let us remark that the definition and properties of the Euclidean Wasserstein distance W2 
discussed above generalizes equally well to any dimension. Tanaka himself j21] showed that 
the Euclidean Wasserstein distance is a non strict contraction for the elastic classical Kac 
model. In the inelastic Kac model, the evolution of the density function / is governed by 
the equation 

^^ QifJ) 



dt 
in which the collision term Q{f, /) is defined by 



(5.1) 



(^,Q(/,/)) = 
for any test function if, where 




2tt 



fiv)f{w) ip{v') -(p{v) 



V COS 6 I cos^l^ — w sin 6 I sin^l^ 



de 

2^ 



■ dvdw 



is the postcoUisional velocity and p > measures the inelasticity. Equation (jS.lj) preserves 
mass but makes the momentum and kinetic energy decrease to at an exponential rate, 



e{f{t))=e-'^'9if) + ie- 



■2f3t 



-2t\ 



<f > with P > given below. In particular, solutions 



to (J5.H) tend to the Dirac mass at 0. 

As in the inelastic Maxwell model discussed above, we start by deriving a contraction 
property for the gain operator Q~^ defined by 

r.27r 



(^,Q+(/,/)) 




dfj 

f{v)f{w) ip{v') -— dv dw 



27r 



Proposition 12. /// and g belong to V2{^), then 



WiiQ'-ifJ),Q^ig,g))< 



r"(|cos^p(P+i) + |sine|2(P+i)) 
Jo 



de 

2^ 



W^{f,g). 
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In terms of solutions f{t) and g{t) to the modified Kac equation (jS-ip with finite initial 
energy only, the above proposition yields, as in previous sections, the bound 

where 

/■27r JQ 

2/? = l- / (|cos^|2(P+i) + |sin^|2(P+i))— >0. 
Jo 27r 

This bound is optimal without further assumptions on the initial data /° and g^ since 
equality holds in the case when < /° >= and g'^ = 60 analogously to previous cases. 

Proof.- Given a vector [v, w) in R^, let Cy^^i denote the curve 

{{v'{e),w'{e)),ee[o,27r]} 

where 

v'{9) = t> cos 6* I cos 6*1^ — w sin 6^ I sin 6*1^ , ^ 

w'{6) = f sin^ I sin6'|P + U7 COS0 I cos^l^. ' 

Let also Ujj,w be the uniform probability distribution on C^^w 

Given V and W two independent random variables distributed according to /, we note 
that Q'^{f,f) is the first marginal on M of E[Z//y_vi/], but also its second marginal by 
symmetry. Then, we have the following result, which is the analogous of Lemmas |21 and 
^2 for this model: 

Lemma 13. Given two vectors {v,w) and {x,y) in M."^, the squared Wasserstein distance 
between the distributions U^^w and Ux,y is bounded by 

(1-2/?) {\v - x\'^ + \w ~ y\^) . 
Proof.- One can transport the curve C^^w onto C^.y by the linear map 

(a, h) 1-^ Tia, b) = —(a costu — b sinu;, a sinu + b cos uj) 
r 



where r = y/v^ + w^, r' = ^J x^ + y^ and uj is the angle between the vectors (f , w) and 
(x, y) in case they do not vanish. We leave the reader discuss the case when either (x, y) 
or (f , w) are zero. Then, analogously to the proof of Lemma one can define a transport 
plan associated to the transport map T to get 

Wl{U,^^M.,^< f \T{a,b)-{a,b)?dU,,^{a,b). 

Furthermore, for all {a,b) in M^, 

|T(a, 6) — (a, 6)p = — (a cos cj — 6 sin cij) — a + — (a sin cj + 6 cos u;) — 6 

= ((7)'-27Cos.; + l)(a2 + 6^) 

If - xP + \w -y\'^ , 2 i2\ 
= 2^2 a +fe • 
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Hence, we deduce 

I 12 I I \2 p 

V — X ~n W — 11 \ I 



y2 _|_ ^2 y^ 27r 

But 

v'{ef + «;'(^)2 = (I cos^|2(P+^) + I sin^|2(P+i))(t;2 + w^) 

by dmi), so that 

r /"St 

L"'o 
which is the bound given by the lemma. D 



/"(|cos^|2(P+^) + |sin^|2(P+i)) — 
Jo 27r 



1^; — x\^ + \w ~ y|^) 



We now continue the proof of Proposition^^ First of all, let (V, X) and (VF, Y) be two 
independent couples of random variables, with V and X distributed according to f , W 
and Y according to (?, optimal in the sense that 

iy|(/, g) =W.[\V- W\^] = E[|X - Fp] . 

Then, by convexity of the squared Wasserstein distance again, it follows from Lemma El 
that 

Wi (E [Uv,w] , E [Uxx] ) < E [Wi{Uv,w. Uxx)] 

< (l-2/3)(E[|F-l^p] +E[|X-r|2]) 

= 2(l-2/5)iy|(/,^). (5.3) 

Next, we remark that the measure Uy^w on "^ has first and second marginals equal 
by symmetry of the curve Cy^w by a 7r/2 rotation. This implies that the first and second 
marginals of E^Wy^^y] on K^ are equal to Q^{f, /), and likewise for the measure E[Z//x,y] 
with marginals Q~^{g,g)- We shall conclude the argument of Proposition ^J by using the 
following general result: 

Lemma 14. // the Borel probability measures /x'- on M are the successive one- dimensional 
marginals of the measure /i* on M^, for i = 1,2 and j = 1, . . . ,N, then 

N 

Proof.- Let vr be a measure on M^ x M^ with marginals /i^ and /i^, optimal in the sense 
that 



Wi{iJ.\ii^)= 11 \v-w\^dn{v,w). 



20 FRANgOIS BOLLEY AND JOSE A. CARRILLO 

Then its marginal ttj on M„, x M^, has itself marginals fi] and /i^, so 



N 



W2{fij,fij)< / / I Vj - Wj I diTj {vj ,Wj). 



The lemma follows by noting that 2. \vj — Wjl"^ = \v — w]"^. D 

In our particular case, Lemma IT^ ensures that 

2W^{Q-'ifJ),Q-'ig,g))<W^{E[Uv,w],^Ux,Y]) 
which concludes the proof of Proposition [T21 taking ()5.3p into account. D 

Appendix: Uniform in time Propagation of fourth order Moments 

In this appendix we derive a uniform propagation of fourth order moments / |f |'^ (^(r, v) dv 
of solutions g to 

^ + V-{gv) = EQ{g,g) (5.4) 

OT 

Q 

where the operator Q(g, g) is defined as in (11.21) for < e < 1 and E = -• 

1 — e"' 
This result has been used in the proof of Theorem [7| 

Proposition 15. If g^ is a Borel probability measure on M.^ such that 

\v\^ g^{v) dv < oo, 
then the solution g to (j5.4j) with initial datum g^ verifies 

sup / \v\'^ g{T,v)dv < oo. 

Proof.- Without loss of generality we can assume that g^, and hence g{T) for all r > 0, 
has zero mean velocity. We let 

m4(r) = / \v\'^ g{T.,v)dv 

denote the fourth order moment of g{T). Then, using the weak formulation of the inelastic 
Boltzmann equation, we have: 

^f^= f V{\v\')-vg{T,v)dv + E [ \v\'Qig{r),g{r)){v)dv. (5.5) 

While the first term in the right hand side is simply 4r7i4(r), the second term is computed 
by 
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Lemma 16. There exist some constants fii and fi2, depending only on e, such that 
\'v\'^Qi9i 9)i''j) dv = — X / \v\'^ g{v) dv + fill / Iv]"^ g{v) dv) 

+ yU2 / / {v ■ w)'^ g{v) g{w) dvdw 



2 



for any probability measure g on M.^ with finite moment of order 4 and zero mean velocity, 

where 

1 1 — e 

A = -(l + 4£-7e^ + 4£^-2e^) and e= 

3^ ' 2 

With this lemma in hand, ()5.5p reads 

dm At) 



(4:-EX]m4{T) + m{T) (5.6) 



dr 

where m{T) is a combination of second order moments, which are bounded in time since the 
kinetic energy is preserved by equation ()5.4|) . Moreover one can check from the expression 
of E and A in terms of e = (1 — e)/2 that 

4-EA = -— ^ -[-l + 2e + e^-4e^ + 2e^] 

36{l — e) 

which is negative for any < e < 1/2, that is, for any < e < 1. By Gronwall's lemma 
this ensures that m^lr) is bounded uniformly in time if initially finite, and concludes the 
argument to Proposition^) D 

Let us remark that identity ()5.6|) is also useful to understand that moments are not 
created by this equation in contrast to the hard-spheres case pTl I22j . In fact, if initially 
moments are infinite, they will remain so. Thus, this is another reason why homogeneous 
cooling states have only certain number of moments bounded (see [Z|). 

We now turn to the proof of Lemma E3 whose result is given in [HI Section 4] only in 
the radial isotropic case, i.e., whenever g{v) depends only on |f |. By symmetry we start 
by writing 



\v\^Q{9,9){v)dv= -— / g{v)g{w)-[\v'\'^+\w'\'^ -Ivl"^ -\w\^] da dvdw 

JlR3 47r Jt^3 7]r3 Js2 z 



where 



, 1, , 1-e, , 1 + e, , 

V = -(v + w)-\ : — [V — w) -\ \v — w\a 

2 4 4 

/ 1/ N 1 ~ e , , 1 + e, , 

w = -[V + w) — [v — w) — \v — w\a. 

Then we introduce the notation 

v+w v—w 1— e , 1+e 

u = 5 U = 5 e = ) e = 1 — e = 

2 2 2 2 

in which 

v' = u + eU + e'\U\a, v' = u - eU - e' \U\(7, v = u + U, w = u-U. 
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Then 

|t;'|2 = \u\'^ + {e^ + e'^)\U\'^ + 2ee'\U\{U-a) + 2e{u-U)+2e\U\{u-a) 
|u,'|2 = \u\^ + {e^ + e'^)\U\^ + 2ee'\U\{U-a)-2e{u-U)-2e'\U\{u-a) 

|t;|2 = |y|2 + |f/|2^2(M-f/) 
|U,|2 = |M|2 + |f/|2_2(M-f/) 

and eventually 

^[\v'\^ + \w'\^ - \v\^ - \w\^] 

= [{e^ + ey - 1] \U\^ + 2(£2 + e'2 _ i)\u\^\u\^ + 4(6^ _ i)(« . [/)2 
+ Ae^e'^ |f/|2 ([/ . a)2 + 4£'2 |f/|2 (tx ■ a)^ 

+ Aee' \U\ [\u\^ + (e^ + e'^)\U\^] (U ■a) + ^ee' \U\ {u ■ U) {u ■ a). 
Integrating with respect to a in 5*2 and taking the identities 

into account, we obtain 

U\v'\^ + \w'\^ - \v\^ - \w\^] ^ = a \U\' + (3 \u\^\U\^ + 7 (m ■ Uf 

g2 2 4:71 



where 



a = (e^ + e'Y - 1 + ^eV^, p = 2[e' + e" - 1 + '^e"] , 7 = 4(e^ - !)• 



Then, by definition of u and f/ in terms of v and w, the identities 

|f/|'^5'(f ) 5'(w) dvdw = -[m^ + ml + 2m2\, 



II 



and 



|-up |f/p (?(f ) g{w) dvdw = -[1114 + ml — 2m2] 



{u ■ UY g{v) g{w) dvdw = -[m^ — m^l] 



hold with 

1714= / \v\'^g{v)dv, m2 = / \v\'^g{v)dv, m^ = I I {v ■ wY g{v) g{w)dvdw 



since g has zero mean velocity. Collecting all terms, we obtain 



1"^! Q{g, g){v) dv = —\m,4 + fiim,2 + fi2'm'2 

where 

A = -lia + /3 + 7) = i(l + 4£ - Te" + 4£3 -2e*), 

o o 
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/ii = -(a + /5-7) and fi^ = -{a - p) 
o 4 

depend only on e, that is, only on e. This concludes the proof of Lemma ITHl D 
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